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ABSTRACT 

/i Orionis was identified by spectroscopic studies as a quadruple star system. Seventeen high preci- 
'^C^'i sion differential astrometry measurements of /i Ori have been collected by the Palomar High-precision 

(•~^ ' Astrometric Search for Exoplanet Systems (PHASES). These show both the motion of the long period 

, binary orbit and short period perturbations superimposed on that caused by each of the components 

■ in the long period system being themselves binaries. The new measurements enable the orientations 

of the long period binary and short period subsystems to be determined. Recent theoretical work 
O , predicts the distribution of relative inclinations between inner and outer orbits of hierarchical systems 

to peak near 40 and 140 degrees. The degree of coplanarity of this complex system is determined, and 
the angle between the planes of the A-B and Aa-Ab orbits is found to be 136.7 ± 8.3 degrees, near 
the predicted distribution peak at 140 degrees; this result is discussed in the context of the handful 
of systems with established mutual inclinations. The system distance and masses for each component 
are obtained from a combined fit of the PHASES astrometry and archival radial velocity observations. 
, The component masses have relative precisions of 5% (component Aa), 15% (Ab), and 1.4% (each of 

^1^' Ba and Bb). The median size of the minor axes of the uncertainty ellipses for the new measurements 

is 20 micro- arcseconds (/ias). Updated orbits for 5 Equulei, k Pegasi, and V819 Herculis are also 
presented. 

^ ' Subject headings: stars:individual(/i Orionis) - stars:individual((5 Equ) - stars:individual(K 

\ Peg) - stars:individual(V819 Her) - binaries: close - binaries:visual - tech- 

niquesiinterferometric - astrometry - stars:distances 



o 

o 



o 



O 



X 



1. INTRODUCTION 



> 

, fi Orionis (61 Ori, HR 2124, HIP 28614, HD 40932) is a quadruple star system that has been extensively studied 
I by radial velocity (RV) and differential astrometry. It is located just North of Betelgeuse, Orion's right shoulde r (left 
P\| on the sky); ^ Ori is a bright star that is visible to the unaided eye even in moderately light-polluted skies. iFrostI 
(|1906D discovered it to be a short period (four and a half day) single-lined spectrosc opic binary; t his was component 
Aa, whose short-period, low mass companion Ab has never been detected directly. |A itken | (119141) discovered it also 
J — ' had a more distant component (B) forming a sub-arcsecond visual binary. Much later, iFekel (1980) found B was itself 
a short-period (4.78 days) double-lined spectroscopic binary, making the system quadruple; these stars are designated 
• • ■ Ba and Bb. Most recently, Fekcl e t al.. (2002) (hereafter F2002) reported the astrometric orbit of the A-B motion, 
double-lined RV orbits for A-B and the Ba-Bb subsystem, and a single-lined RV orbit for the Aa-Ab subsystem. 
F2002 estimate the spectral types as A5V (Aa, an Am star), F5V (Ba), and F5V (Bb), though they note these are 
classifications are less certain than usual due to the complexity of the system. For a more complete discussion of the 
' history of ^ Ori, see F2002. 

Until now, astrometric observations have only been able to characterize the long period A-B motion, lacking the 
precision necessary to me asure the astrometric perturbat ions to this orbit caused by the Aa-Ab and Ba-Bb subsystems. 
The method described bv lLane and Muterspaughl ()2004n for ground-based differential astrometry at the ~ 20 yitas level 
for sub-arcsecond ("speckle") binaries has been used to study ^ Ori during the 2004-2007 observing seasons. These 
measurements represent an improvement in precision of over two orders of magnitude over previous work on this 
system. 

The goal of the current investigation is to report the center-of-light (photocenter) astrometric orbits of the Aa-Ab 
and Ba-Bb subsystems. This enables measurement of the coplanarities of the A-B, Aa-Ab, and Ba-Bb orbits. The 
masses and luminosity ratio of Aa and Ab are measured for the first time. Also presented are updated orbits for the 
PHASES targets 5 Equ, n Peg, and V819 Her. 
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Astrometric measurements were made at the Palomar Testbed Interferometer (PTI; IColavita et al.lll999|) as part o f 
the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) program ( Muterspaugh et al]l2006cf ). 
PTI is located on Palomar Mountain near San Diego, CA. It was developed by the Jet Propulsion Laboratory, California 
Institute of Technology for NASA, as a testbed for interferometric techniques applicable to the Keck Interferometer 
and other missions such as the Space Interferometry Mission (SIM). It operates in the J (1.2^m), H (1.6/im), and K 
(2.2/ini) bands, and combines starlight from two out of three available 40-cm apertures. The apertures form a triangle 
with one 110 and two 87 meter baselines. 

2. OBSERVATIONS AND DATA PROCESSING 
2.1. PHASES Observations 
2.1.1. Instrumental Setup 

/i O ri was observed with PTI on 17 nights in 2004-2007 with the observing mode described in lLane and Muterspauglil 
(|2004f ). Starlight is collected from two apertures, coUimated, and sent to a central beam combining facility. There, 
the light from each telescope reflects from movable mirrors (delay lines) whose position is constantly varied to account 
for sidereal motion and to track atmospheric index of refraction variations. After this first set of delay lines, a beam- 
splitter is used to divide the light from each telescope; ~ 70% of the light is sent to an interferometric beam combiner 
that monitors a single fringe from any star in the field at rapid (10-20 millisecond) time-scales to measure fringe phase 
variations caused by the atmosphere, and provide feedback to the mai n delay lines. This pro cess phase-stabilizes 
the other ^ 30% of the light (a technique known as phase referencing, iLane and Colavitall2003l ). which is sent to a 
second interferometric beam combiner that can add an additional variable delay of order 250 /im to the light from one 
telescope. This variable delay is modulated with a triangle waveform, scanning through interferograms from all stars 
within the subarcsecond field of view. These interferogram scans are the observables used for PHASES astrometry. 

2.1.2. Data Reduction 

Modifications to the data processing algorithm since the original report are given by iMuterspaugh et ahl ()2005f ) 
and have been incorporated in the current study. Interferogram templates are fit to each scan, forming a likelihood 
function of the separations of interferograms formed by components A and B. A grid of differential right ascension 
and declination is formed, and a likelihood surface is mapped onto this grid by converting delay separation to 
differential right ascension and declination. That surface is coadded over all the scans of /i Ori within the night 
(typically ~ 1000 scans or more). The deepest minimum in the surface corresponds to the binary separation, while 
the width of that minimum determines the uncertainty ellipse. Due to the oscillatory nature of the interferograms, 
other local minima can exist; these "sidelobes" are separated by the interferometer's resolution ~ A/B ~ 4 mas {B is 
the separation between the telescopes, and A the wavelength of starlight), an amount much larger than the width of 
an individual minimum. The SNR can be increased by coadding many scans and by earth-rotation synthesis, which 
smears all but the true minimum, and the true minimum then can be established. Only those measurements for which 
no sidelobes appear at the Aa contour of the deepest minimum are used in orbit fitting. 

All measurements have been processed using this new data reduction pipeline. The measurements are listed in Table 
m in the ICRS 2000.0 reference frame. 

2.1.3. Technique Upgrades 

Data since mid-2006 have benefited from the use of an automatic alignment system and longitudinal dispersion 
compensator; the affected data points are noted in Table [TJ These modifications reduced the throughput of the 
astrometry setup; to compensate, a 50 Hz phase tracking rate was sometimes used, whereas observations previous to 
these changes utilized 100 Hz tracking for monitoring the atmosphere. 

Drifts in optical alignment may result in variable pupil sampling at the interferometer apertures, changing the 
effective interferometric baseline. To minimize this potential systematic, a continuous realignment system has been 
developed. A red laser is coaligned with the starlight and reverse propagated through the interferometer. Four-percent 
reflective pellicle beamsplitters are placed near the focuses of the interferometer telescopes to extract this tracer beam 
and redirect to a camera where the pupil is reimaged. The angles of the first flat mirrors receiving incoming starlight 
in the beam combining lab are continually adjusted by a closed-loop feedback system to hold the laser spot on the 
camera at the telescope. 

The path compensation for the geometric delay at PTI has been done with movable mirrors in air, which has a 
wavelength-dependent index of refraction. The fringe packets of astrophysical sources are dispersed by an amount that 
depends on the difference in air paths between arms of the interferometer; this changes the shape and overall location 
of the fringe packets. To compensate, two prisms are introduced in each of the interferometer's two arms. The set in 
one arm is static. The other pair are slid relative to each other along the slope of the prism to introduce a variable 
amount of glass dispersion whose shape is opposite that of air to high order. This flattens the variability of delay versus 
wavelength. The system is calibrated by setting the telescope siderostats into a retroreflecting mode, using an internal 
white-light source to form interferograms, which are detected with a low-resolution spectrometer (5 elements across K 
band), and measuring the offsets between the interferograms as a function of prism location. During observations, the 
prism position is continuously changed with an open loop control calculated from the locations of the delay lines. 

There is insufhcient new data to establish the degree two which these instrumental changes might be reducing excess 
data scatter or to establish a relative weighting between data subsets. No large discontinuities in the orbital motions 
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TABLE 1 

PHASES DATA FOR ^l OrI 
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Note. — All quantities are 


in the ICRS 2000.0 reference frame. 


The uncertainty values presented in 


these data have 


all been 



scaled by a factor of 1.73 over the formal (internal) uncertainties within each given night. Column 1 is the heliocentric modified Julian 
date. Columns 2 and 3 is the differential right ascension and declination between A and B, in milli-arcseconds. Columns 4 and 5 
arc the la uncertainties in the minor and major axis of the measurement uncertainty ellipse, in micro-arcseconds. Column 6, 0c, is 
the angle between the major axis of the uncertainty ellipse and the right ascension axis, measured from increasing differential right 
ascension through increasing differential declination (the position angle of the uncertainty ellipse's orientation is 90 — 0c). Columns 7 
and 8 are the projected uncertainties in the right ascension and declination axis, in micro-arcseconds, while column 9 is the covariance 
between these. Column 10 is the number of scans taken during a given night. Column 11 is 1 if the longitudinal dispersion compensator 
was in use, otherwise. Column 12 is 1 if the autoaligncr was in use, otherwise. Column 13 represents the tracking frequency of 
the phase-referencing camera. The quadrant was chosen such that the larger fringe contrast is designated the primary (contrast is a 
combination of source luminosity and interferometric visibility). 

are seen between pre- and post-upgrade subsets for /i Ori or the other PHASES targets, suggestmg the subsets are 
compatible. For the purposes of the current investigation, the PHASES observations are treated as a single data set 
with equal weighting on observations from before and after these upgrades. 

2.1.4. The PHASES Astrometnc Orbits 

The differential astromctry measurements are listed in Table[l] in the ICRS 2000.0 reference frame. The existence of 
data scatter beyond the level estimated by formal uncertainties from the PHASES analysis algorithm was determined 
by model fitting the PHASES data alone. Model fitting was performed with standard sum of squared residuals 
minimization, slightly complicated by the two dimensional nature of the uncertainty ellipses but still straight-forward 
to carry out. The limited number and time span of the PHASES observations prevents an independent fitting of that 
data set to a 4-body, 3-Keplerian model to determine potential noise excess. Thus the A-B period and Aa-Ab and 
Ba-Bb periods, eccentricities, and angles of periastron passages were fixed at the values reported in F2002 (in the 
case of Ba-Bb, which had zero eccentricity in F2002, the angle of periastron passage is undefined and fixed at zero). 
The minimum does not equal the number of degrees of freedom. An excess noise factor of 1.73 is found, and the 
PHASES uncertainties reported in Table [T] have been increased by this amount over the formal uncertainties. The 
rescaled (raw) median minor- and major-axis uncertainties are 20 (11) and 347 (200) /ias. The rescaled (raw) mean 
minor- and major-axis uncertainties are 27 (16) and 668 (386) /ias. 

2.2. Previous Measurements 

Previous differential astrometry measurements of /i Ori are tabulated in Table 5 of F2002. These have been included 
in the current fit, with identical weightings as assigned by that investigation, though it is noted that the text contains 
a typographical error, and the p unit uncertainty ap should be 0.024, rather than 0.0024 mas. The time span of 
these measurements is much longer than that of the PHASES program and aids in solving the long period A-B orbit, 
which also lifts potential fit parameter degeneracies between that orbit and those of the short period subsystems. 
Measurements marked as 3cr outliers by that investigation are omitted, resulting in 80 measurements each of separation 
and position angle being used for fitting. Ten new measurements have been published since that investigation and 
are listed in Table [2] with weights computed with the same formula as used in F2002. Two of these measurements are 
found to be 3(T outliers. In total, there are 88 measurements of separation and position angle used in fitting. 

F2002 also present radial velocity observations of components Aa, Ba, and Bb. Those measurements are included in 
the present fit, with weightings as reported in Tables 2, 3, and 4 of that paper. Measurements marked as 3a outliers 
by that investigation have not been included in the present analysis. In total, 442 velocities — 220 for Aa and 111 for 
each of Ba and Bb — are used in fitting. 



3. ORBITAL MODELS 
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TABLE 2 

New Non-PHASES Astrometry for ^l Ori 



Besselian Year p 9 Weight Outlier Reference 

1991.8101 0.330 31.8 0.1 1 Fabricius et al, 

1996.8986 0.306 13.4 8.7 Horch et al, 

1999.0153 0.200 14.8 10.1 1 Horch et al. 



1999.0153 0.196 13.3 9.7 1 Horc 



1 et a 



1999.0153 0.203 14.9 10.4 1 Horch et a 



1999.8915 0.150 11.7 5.8 1 Horch_et_al. 

1999.8915 0.145 12.5 5.5 1 Horch et al. 

1999.8915 0.154 11.3 6.1 1 Horch et al. 

2000.7653 0.081 359.4 2.0 Horch et al. 



2005.1331 0.179 26.4 0.2 1 Scardia et al. 




Note. — The 10 new astrometry measurements published since 
F2002 for i^l Ori. Column 1 is the epoch of observation in years, column 
2 is the A-B separation in arcscconds, column 3 is the position angle 
east of north in degrees, and column 4 is the measurement weight on 
the same scale as F2002. ((5RA — ps\n6, (5Dec — pcosQ.) Column 5 
is if the measurement is a 3(T outlier not used in fitting, 1 otherwise. 
Column 6 is the original work where the measurement was published. 

The apparent motions of the centers-of-hght (photocenters) of the A = Aa — Ab and B = Ba — Bb subsystems rela- 
tive to each other are given by the model 



2/obs = f A-B 

^ MAb/MAa - ^Ab/^Aa 

(1 + MAb/MAa) (1 + iAb/^Aa) 
Msb/MBa — L^hl Lb'a 



^Aa-Ab 



(1 + MBb/MBa) (1 + iBb/^Ba) '"^^ 

corresponding to a four-body hierarchical dynamical system (HDS). The quantities M are component masses and 
L are component luminosities, and each of the summed vectors is determined by a 2-body Keplerian model. This 
model is used to fit the astrometric data; note that the total masses from the Aa-Ab and Ba-Bb orbits also show 
up as component masses for the A-B orbit, linking them, and that the mass ratios and luminosity ratios appear as 
additional parameters, degenerate with each other. The radial velocities are fit by a simple superposition of individual 
Keplerians; these determine mass ratios, and the luminosity ratios become nondegenerate parameters. The luminosity 
ratios are primarily constrained by the K-band PHASES data; the other astrometric data are not precise enough to 
detect the subsystem motions. 

The combined simultaneous fit to all data sets has 26 free parameters and 626 degrees of freedom. The parameters 
used are listed in the top half of Table [3] with their associated best-fit values and Icr uncertainties. It should be noted 
that while the Ba-Bb eccentricity has been allowed to vary as a free parameter, the best fit value is consistent with 
zero and could have been fixed; the other parameters are not changed significantly by doing so. Quantities of interest 
derived from those parameters are listed in the second half of that table, with corresponding uncertainties derived 
from first-order uncertainty propagation. The apparent center-of-light wobbles of the Aa-Ab and Ba-Bb subsystems 
are plotted in Figure (TJ the A-B and RV orbits were plotted in F2002 and are not significantly different in the present 
model. 

Including PHASES measurements in the fit introduces the ability to evaluate the inclination and luminosity ratio of 
the Aa-Ab system and angles of the nodes of the Aa-Ab and Ba-Bb systems, quantities that were entirely unconstrained 
in the F2002 study. The A-B angular parameters have much smaller uncertainties than in F2002 (the ^ab, *ab, and 
wab uncertainties are reduced by 13 x, 6x, and 5x, respectively). Uncertainties in the A-B period, eccentricity, and 
epoch of periastron passage are improved by a factor of 2 or more. Most other fit parameters are constrained only 
marginally better than in F2002. 

3.1. Relative Orbital Inclinations 

In systems with three or more stars, studying the syst em coplanarity is of interest for understanding the formation 
and evolution of multiples (jSterzik and Tokovinin|[2"002f ). To determine this without ambiguity, one must have both 
visual an d RV orbital solu tions for pairs of interest. Previously, this was available only in six triples (for a listing, see 
[Muterspai lgh et al.|[2006al ). The reason mutual inclination measurements have been rare is due to the observational 
challenges these systems present — RV signals are largest for compact pairs of stars, whereas astrometry prefers wider 
pairs. The "wide" pair must be studied with RV and thus have an orbital period (and corresponding separation) as 
short as the two-component binaries that are already challenging to visual studies. The "narrow" pair is even smaller. 
The mutual inclination between two orbits is given by: 

cos <& = coszi cosi2 -I- sinii sin ^2 cos (fii — ^2) (2) 

where ii and 12 are the orbital inclinations and Oi and ^2 are the longitudes of the ascending nodes. If only velocities 
are available for one system, the orientation of that orbit is unknown (even if it is eclipsing, the longitude of the 
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TABLE 3 
Orbital parameters for ^l Ori 





Lxa. > LAb, 


LAa > i-Abi 


LAa < i-Abi 


^Aa < i-Ab) 


Parameter 


Lsa. > ^Bb 


-i'Ba < ^Bb 


-^'Ba > -^'Bb 


^Ba < ^Bb 



■Pab (days) 
CAB 

*AB (degrees) 
UJAB (degrees) 
Tab (MHJD) 
f^AB (degrees) 
-PAaAb (days) 
EAaAb 

lAaAb (degrees) 
t^AaAb (degrees) 
TAaAb (MHJD) 
J^AaAb (degrees) 
-PBaBb (days) 

fiBaBb 

JBaBb (degrees) 
t^BaBb (degrees) 
TeaBb (MHJD) 
f^BaBb (degrees) 
Ma (Mq) 
Mb (Mo) 

MAb/MAa. 
-MBb/Mea 
I^Ab/ Las. 
^Bb/i-Ba 
d (parsecs) 
Vq (kms-1) 



6813.8 
0.7410 

96.028 
36.712 
46090.7 
204.877 
4.4475849 
0.0037 
47.1 
304 

43739.69 
50.5 

4.7835349 
0.0016 
110.71 
217 

43746.40 
111.3 
3.030 

2.746 

0.274 

0.9764 



0.765 
46.11 

42.548 



723.86 
±1.2 
±0.0011 
±0.028 
±0.066 
±1.0 
±0.011 
±1.2 X 10 
±0.0014 
±9.0 
±21 
±0.26 
±3.7 
±3.0 X 
±0.0014 
±0.73 
±47 
±0.63 
±3.9 
±0.069 
±0.038 
±0.051 
±0.0022 
±0.040 
±0.055 
±0.28 
±0.027 



723.86 



723.58 



723.58 



10" 



50.0 



±8.1 



231.7 ±3.8 



291.3 ±3.9 



A^Aa (Mq) 2M IOTT 

MAh (Mq) 0.652 ±0.097 

Mfla (Mq) 1.389 ±0.019 

A^Bb (Mq) 1.356 ±0.019 

*AB-AaAb (degrees) 136.7 ±8.3 

*AB-BaBb (degrees) 91.2 ±3.6 

*AaAb-BaBb (degrees) 84.6 ±4.9 

OAB (mas) 273.7 ±2.1 

aAB (AU) 12.620 ±0.057 

aAaAb,coL (^as) 358 ±84 

OAaAb (mas) 1.661 ±0.016 

lAaAb (AU) 0.07659 ±0.00058 

OBaBb.cOL (/ias) 102 ±30 

OBaBb (mas) 1.688 ±0.013 

OBaBb (AU) 0.07780 ±0.00036 

TT (mas) 21.69 ±0.13 

ifAa 1-03 ±0.26 
KAb > 4.58 

Kbs. 1.72 ±0.26 

A'Bb 2.02 ±0.26 

L^.Aa 8.3 ±2.0 

Lk, Ab ±0.33 

Lk,b>, 4.4 ±1.1 

Lji Bb 3.35 ±0.82 



246 ±0 



84.5 ±3 
125.1 ±6 



1.99 
1.75 



3.45 
4.3 



111.3 ±4.0 

0.259 ±0.039 

0.738 ±0.061 

089 0.773 ±0.055 



291.3 ±4.0 



233 ±0 



2.408 ±0.092 
0.623 ±0.075 



52.2 ±6.1 
91.2 ±3.8 
126.2 ±5.9 



84.5 
82.2 



±3 
±4 



NiiTK. — Orbital parameter 
the last reported digit fron 
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±30 










1.64 


±0.26 










1.96 


±0.27 






±0 


26 


1.73 


±0.26 


1.98 


±0.26 


±0 


26 


2.01 


±0.26 


1.75 


±0.26 






4.8 


±1.2 










3.52 


±0.86 






±0 
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4.4 


±1.1 


3.47 


±0.84 


±1 





3.38 


±0.82 


4.3 


±1.0 


cllipB 


□ s ind 




ametcr that 




by less tha 



088 



neter uncertainties have been increased by a factor of Xr ' 



1.08 



; for p. Ori. In the second, third, and fourth solutic 
the previous model. In the combined fits, all par 

(though the of the combined fit is artificial due to rescaling the uncertainties of the individual data sets, this reflects the degree to which the 
data sets agree with each other). The first solution is strongly preferred as it produces masses and luminosities that are correlated; the second is 
also possible because the stars Ba and Bb are very similar. The third and fourth solutions require an unlikely luminosity for component Ab, given 
its mass, and are not preferred. a^oL *s the semimajor axis of the motion of the center-of-light of a subsystem, at K-band. L-^/Ly is the K-band 
luminosity ratio between components X and Y. K is the K-band absolute magnitude, L is K-band luminosity, in solar units. For the first two 
solutions, the best fit solution yields H i^ infinite; a lower limit is determined by setting ^Ab/^Aa ^'-^ the upper limit of its Ict uncertainty range, 

ascending node is unknown). If velocities are unavailable for a given orbit, there is a degeneracy in which node is 
ascending — two values separated by 180 degrees are possible. This gives two degenerate solutions for the mutual 
inclination (not necessarily separated by 180 degrees). 

Even when a center-of-light astrometric orbit is available for the narrow pair, there can be a degeneracy between 
which node is ascending and the luminosity ratio. Having found one possible luminosity ratio LAb/^Aa — Li, it can 
be shown that the other possible solution, corresponding to changing the ascending node by 180 degrees, is given by 

2R + RLi—Li . . 

= l + 2L,-R 



where R is the mass ratio MAb/AfAa- In previous studies, support data has been available to lift that degeneracy. For 
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Aa-Ab Center of Light Motion Ba-Bb Center of Light Motion 

Z I 1 1 1 r 




time (days) time (days) 

Fig. 1. — The astrometric orbits of /i Ori Aa-Ab and Ba-Bb, phase- wrapped about their respective orbital periods. Phase zero is at the 
epoch of periastron passage (T), and each plot is repeated for two cycles (and each measurement is plotted twice) to allow for continuity at 
all parts of the graph. In both cases, the motions of the A-B system and of the other subsystem have been removed. The projection axis 
shown for each is 159 degrees East of North (equivalent to position angle 291 degrees), well aligned with the minor axis of many PHASES 
observations. The plotted uncertainties are those projected along this axis, and have been increased by a factor of 1.73 over the formal 
uncertainties, to reflect excess noise within the PHASES set. On the left is the motion of Aa-Ab; for clarity, only those observations with 
rescaled and projected uncertainties less than 200 fias have been plotted; on the right is the motion of Ba-I3b, for which the cutoff was at 
30 /jas. 

example, in the V819 Her system (jMuterspaugh et aLll2006al ). the two possible luminosity ratios were 0.26 and 1.89. 
The eclipsing nature of the Ba-Bb pair lifted that degener acy becau se it helped establish that the luminosity ratio is 
much less than 1. Similarly, in the k Peg system ( Muterspaugh et al. 2006b), the spectra used for RV also show that 
component Bb is much fainter than Ba, again lifting the degeneracy (a luminosity ratio either nearly zero or 1.9 were 
both possible). 

For a quadruple system, there are as many as four degenerate fit solutions. For /i Ori, a global minimum is found 
with longitude of the ascending node f^AaAb — 231.7 ± 3.8 degrees and LAb/^Aa = 0.738 ± 0.061. The alternative 
pair of these parameters solving eq. [3] would force the luminosity ratio to a negative value (but close to zero within 
uncertainties). However, searching for solutions with nonnegative luminosity ratios near zero yields a fit solution with 
only shghtly larger x^, and the node at roughly 180 degrees difference (f^AaAb = 50.5±3.7 and LAh/LAa = 0±0.040). 
All parameters other than the node angle and luminosity ratio vary between the two models by amounts less than 
the fit uncertainties. The two solutions find MAbZ-^Aa = 0.259 ± 0.039 or 0.274 ± 0.051 respectively; given the rough 
scahng L oc (jSrnith 1983), it is very likely that the larger luminosity ratio is incorrect. The larger luminosity ratio 
would also imply Ab is as bright as Ba or Bb. However, component Ab is not observed in the spectrum while Ba and 
Bb are. This is suggestive that Ab is faint, though it is also possible to explain this lack of Ab lines by postulating it 
is rapidly rotating. However, given that the other stars are not rapid rotators, there is little evidence to support Ab as 
a rapid rotator. It is concluded that the luminosity ratio near zero is strongly preferred, despite the slightly worse 
fit. Both fits are reported in Table [3l but the rest of the discussion in this paper refers only to the preferred solution 
for Aa-Ab. This degeneracy can be fully lifted by a single epoch image with a closur e phase capable interfe rometer 
with sufficient angular resolution (such as the Navy Prototype Optical Interferometer. [Armstrong et al.iri998f) . 

For each of the Aa-Ab solutions, there exists two solutions for the Ba-Bb pair. In these cases, no negative luminosity 
ratios are found; the degeneracy is perfect and the of the fits are identical. f^BaBb differs by 180 degrees in the two 
orbits, and the luminosity ratio LBb/^Ba switches between being larger (at f^BaBb = 291 degrees) or smaller (at 111 
degrees) than unity; all other parameters remain unchanged. Both solutions are near unity, and the stars have very 
similar masses (MBb/-^Ba = 0.9764 ± 0.0022). While the solution for which Bb is less luminous than Ba is slightly 
more consistent because it correlates to the mass ratio, it is conceivable that the other solution is correct. Thus, the 
solution for which Bb is less luminous than Ba is slightly preferred, but not as conclusively as for the Aa-Ab case, 
where the differences between the stars are more significant. Thus, both possibilities are considered in the remainder 
of this paper. 

3.2. Evidence for Kozai Cycles with Tidal Friction? 

Of p articular int erest is the potential for Kozai oscillations between orbital inclination and eccentricity in the narrow 
pairs (|Kozai|[l962f ). which can affect the orbital evolution of the system. These occur independent of distances or 
component masses, with the only requirement being that the mutual inclination is between 39.2 and 180 — 39.2 = 140.8 
degrees. Other effects that cause precession can increase the value of this critical angle. 
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iFabrvckv and Tremaind (|2007f ) predict a buildup of mutual inclinations near 40 and 140 degrees by the combined 
effects of Kozai Cycles with Tidal Friction (KCTF) for triples whose short-period subsystems have periods between 3 
and 10 days. The mutual inclination of /i Ori AB-AaAb is near 140 degrees, leading one to consider if this trend is 
starting to be seen. The systems with unambiguous mutual inclinations break down as follows: 

• Five systems are outside of the 3-10 day inner period range. These systems do not meet the criteria to 
be included in testing the buildup prediction: 

1. V819 Her ($ = 26.3 ± 1.5 degrees, 2.23 d; see g]), 

2. Algol = 98.8 ± 4.9 degrees, 2.9 d; iLestrade et aIlll993HPan et al.lll993[ ). 

3. ?7 Vir ($ = 30.8 ± 1.3 degrees, 72 d; iHummel et al|[200l . 

4. C Uma ABC ($ = 132.1 degrees, 670 d; lHeintzl[l996| ). and 

5. e Hya ABC ($ = 39.4 degrees, 5500 d; Eihitz 1996). 

These fall outside the 3-1 day range of inner-binary periods applicable to the prediction in 
IFabrvckv and Tremaind (|2007f ). However, it is worth noting the mutual inclinations of ^ UMa and e Hya are 
near 140 and 40 degrees respectively and in V819 Her and rj Vir the values are outside the 40 — 140 degrees 
range, so neither would have been predicted to undergo Kozai cycles or KCTF. Algol has a nearly perpendicular 
alignment, though the dynamics of Algol are different due to quadrupolc distortions in the semidetached stars; 
Algol's alignment has been explained by Egglcton and Kiscleva-Eggleton (2001). 

• No systems are outside the 40-140 degrees range, while also in the 3-10 day inner period range. 

• Two systems are between 40 and 140 degrees and in the 3-10 day inner period range. These systems 
would not support the KCTF-driven buildup near 40 and 140 degrees: 

1. fi Ori AB-BaBb ($ = 91.2 ± 3.6 or 84.5 ± 3.6 degrees are possible, 4.78 d) and 

2. 88 Tau AaAb-AblAb2 ($ = 82.0 ± 3.3 or 58 ± 3.3 degrees are possible, 7.89 d; iLane et al.lf2007l ). 

While mutual inclination degeneracies continue to exist in both systems, all possible values are in this range. 

• Three systems are near 40 or 140 degrees, and in the 3-10 day inner period range. These systems 
would appear to support the KCTF prediction: 

1. n Ori AB-AaAb = 136.7 ± 8.3 degrees, 4.45 d), 

2. 88 Tau AaAb-AalAa2 ($ = 143.3 ± 2.5 degrees, 3.57 d; iLane et al.|[2007[) . and 

3. K Peg ($ = 43.4 ± 3.9 degrees, 5.97 d; see g]). 

In total, 3 of the 5 systems meeting the criteria for testi ng the KCTF prediction do app ear near the peak points of 
40 and 140 degrees. Following equations 1, 22, and 35 in IFabrvckv and Tremaing (|2007| ). in the presence of general 
relativity (GR) precession one expects the critical angles for fj, Ori AB-AaAb and k Peg to be increased from 39.2 
degrees to ~ 68 {tujgr = 2.3) and ^ 54 degrees {tlugb. — 1-3), respectively, while for 88 Tau AaAb-AalAa2 GR 
precession dominates no matter the inclination {tujgr = 17). Thus, Kozai oscillations are suppressed by precession 
in these systems' current states. However, it is possible these were present at earlier stages in the systems' histories 
and their current configurations were still reached through KCTF — the inner binaries may have originally been more 
widely separated, in which case GR effects would have been reduced. 

Alternatively, n Ori AB-BaBb and 88 Tau AaAb-AblAb2 both lie well within the range of predicted Kozai cycles, 
even including GR precession (which raises the critical angles to ~ 58 {tujgr = 1-6) and ~ 50 degrees {tlogr — 0.9), 
respectively) . 

Several more systeni s with double visual orbits but lacking RV for at least one subsyster n component are listed by 
ISterzik and TokovininI (|2002D . Two more (HP 150680 and HD 214608) are mentioned bv lOrlov and Petroval ()2000l ) 
and one more (HD 108500) bv lOrlov and Zhuchkovl ()2005[ ). Of these, the inner pair in HD 1506 80 is doubtful and 
listing HD 214608 as havi ng a dou ble visual orbit appears to be in error. In the paper cited bv lOrlov and Petroval 
()2000l ) for HD 214608, Du auennovl (l987) reveal it to be a double spectroscopic system, but point out that the the 
visual elements of the inner pair are unconstrained. The value of 150 degrees for the node seems to have been taken 
as nominal from the outer system (whose true ascending node is 180 degrees different). Because the nodes cannot be 
distinguished in the visual-only pairs, two values of the mutual inclinations are equally possible for each. Furthermore, 
all have inner systems with periods longer than 300 days. Thus, these cannot provide further direction on testing the 
KCTF prediction. 
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Fig. 2. — Stellar evolution models predicting isochrones for the elements of the ^ Ori system show disagreement between component Aa 
with Ba and Bb. Curve labels give values of log(Age/years). Component Aa provides the most leverage for determining the system age, 
though its evolution may have been altered by tidal friction. 

3.3. Masses and Distance 

Components Ba and Bb are each determined to 1.4%, Aa to 5%, while the lowest mass member, Ab, is uncertain 
at the 15% level. The individual masses of Aa and Ab have not been previously determined; this study enables the 
exploration of the natures of those stars. Both are members of star classes of interest: Aa is of spectral type Am, and 
Ab has a mass in the range of late K dwarfs. Through their physical association, it can be assumed both are co-evolved 
with Ba and Bb, each within the mass range for which stellar models have been well calibrated through observation. 

The masses of Ba and Bb are determined only slightly better (less than a factor of 2 improvement) over the previous 
study by F2002. Similarly, the distance is determined t o 0.6% , a slight improvement. The Hipparcos f Perrvman et al] 
[1997) based parallax values discussed by Sodcrhielm (119990 (21.5 ± 0.8 mas in the original evaluation, revised to 
20.8 ± 0.9 mas when the binary nature was considered) are consistent with, but less well constrained, than the current 
value of 21.69 ± 0.13 mas. 



3.4. Component Luminosities 

The 2MASS K-band magnitude for fj. Ori is rnxotai = 3.637±0.260 (jSkrutskie et al.ll2006f) . F2002 gives the difference 
between the luminosities of A and B in several bands. Unfortunately, none of these were taken near the K band (2.2 
^m) where PTI operates. However, a Keck adaptive optics image of fj, Ori was obtained on MJD 53227 with a narrow 
band H2 2-1 filter centered at 2.2622 microns. The A-B differential magnitude in this band is toa — rns — Am^s = 
—0.073 ± 0.007 magnitudes; this measurement is reported for the first time here. 

The combined orbital fit provides the system distance d = 46.11 ± 0.28 parsecs and the luminosity ratios. The 
combined set of rnxotai, Attiab, d, iAb/^Aa, and LBb/iBa determines the component luminosities. Using first order 
error propagation, the K-band luminosities are Lk, Aa = 8.3 ±2.0 solar and less than a third solar for Ab. Components 
Ba and Bb have K-band luminosities of cither Lj^ Ba = 4.4 ± 1.1 and L^, Bb = 3-35 ± 0.82 or Ba = 3.45 ± 0.84 and 
Lk. Bh = 4.3 ± 1.0 solar. Absolute magnitudes are also given in Table[31 In each case, the uncertainty in the apparent 
magnitude WTotai dominates. 

3.5. System Age and Evolutionary Tracks 

The masses and abs olute K-band magn itudes for the components in this system can be compared to the stellar 
evolution models from Girardi et al.l (|2002D . As in F2002, an abundance of Z = 0.02 is assumed. Figure [5] shows the 
mass vs. K-band magnitudes for several isochrones downloaded from http:/ /pleiadi. pd.astro.it/ and the values derived 
for the components of /i Ori. As the most massive component, the properties of Aa provide the strongest constraints 
on age, being most consistent with isochrones in the age range of 10^ — 10*'^ years. This is not entirely consistent with 
the properties of Ba and Bb, though close. Of course, if KCTF has played a significant role in the orbital evolution 
of this system, one wonders whether stellar evolution models for single stars are really applicable to these stars. One 
would anticipate the evolution of Aa as being affected by tidal forces because it is part of the subsystem near the 
predicted 140 degree "pile-up". Thus, one would rely more on Ba and Bb for system age determination, indicating an 
age over 10^ years. 

4. UPDATED ORBITS 

During the course of this inv estigation, a sign error wa s found in the analysis software that was used to compute 
orbital solutions for S Equulei (|Muterspaugh et al.l 120051 ). n Pegasi ([Muterspaugh et al.ll2006"bl ). and V819 Herculis 
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(|Muterspaugh et al.ll2006al ) . This error affected fits to the radial velocity data only, with the result that the descending 
node was misidentified as the ascending, and the angle of periastron passage is off by 180 degrees. Because the software 
was self-consistent, this has no impact on the mutual inclinations derived. Additionally, the finite travel time of light 
across the wide orbit is included in the analysis, and this amount was thus incorrect by the same sign error, but the 
light travel time correction has only a small impact on those models. Both errors have since been corrected. 

Twenty-three new observations of 5 Equ and 11 of k Peg have been collected since those initial investigations and 
are presented in Table [D Also presented in Table 2] is the complete set of 34 V819 Her PHASES observations with 
the 10 measurements taken during eclipses marked; the previous investigation used less precise methods for predicting 
eclipse times, so the flagged measurements have changed. Measurements made during eclipse are not used in fitting. 

The new analysis makes use of the V819 Her Ba-Bb inclination constraint derived by eclipse lightcurves, a feature 
not included in the previous study. When computing x^, an additional term (isaBb - *BaBb,eciipso)^/o'f_BaBb, eclipse is 
added to the sum, where iBaBb. eeiipse — 80.63 a nd aj RaBb, peiipse = 0.33 d egrees are the value and uncertainty of the 
Ba-Bb inclination from the lightcurve studies of Ivan Hamme et al.l ()1994( ). Note that this measurement results from 
an entirely independent data set. This added constraint lessens covariances between orbital elements. 

The corrected and updated orbital solutions are presented in Tables [5] and [6l which are fit to the complete set of 
PHASES observations and the other astrometric and RV observations tabulated in those previous papers. (A few 
new measurement from speckle interferometry are available for each system since those investigations. These have 
little impact on the orbital solutions are not included in the present fit to avoid over-complicating this update.) The 
updated reweighting factors for the PHASES uncertainties for each data set are 3.91 for S Equ, 7.93 for k Peg, and 
2.0 for V819 Her. Alternatively, the noise floor for k Peg is now found at 161 /ias. 

5. CONCLUSIONS 

The center-of-light astrometric motions of the Aa-Ab and Ba-Bb subsystems in /x Ori have been constrained by 
PHASES observations. While four degenerate orbital solutions exist, two of these can be excluded with high reliability 
based on mass-luminosity arguments, and the fact that Ab is not observed in the spectra. Ba and Bb are stars of a 
class (mid-F dwarfs) whose properties have been well established by studying other binaries. Their association with 
Aa and Ab, which are members of more poorly studied classes (Am and late K dwarfs) allows a better understanding 
of those objects in a system which can be assumed to be coevolved. The orbital solution finds masses and luminosities 
for all four components, the basic properties necessary in studying their natures. 

Complex dynamics must occur in /i Ori. The Ba-Bb orbital plane is nearly perpendicular to that of the A-B motion, 
and certainly undergoes Kozai-type inclination-eccentricity oscillations. It is possible that the mutual inclination of 
the A-B pair and Aa-Ab subsystem is a result of KCTF effects over the system's evolution. 

Finally, it is noted that the orbits in the fj, Ori system are quite non-coplanar. This is in striking contrast with the 
planets of the solar system, but follows the trend seen in triple star systems. With the solar system being the only one 
whose coplanarity has been evaluated, it is difficult to draw conclusions about the configurations of planetary systems 
in general. It is important that future investigations evaluate the coplanarities of extrasolar planetary systems to 
establish a distribution. Whether that distribution will be the same or different than that of their stellar counterparts 
may point to similarities or differences in star and planet formation, and provide a key constraint on modeling multiple 
star and planet formation. 
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TABLE 4 

New phases data for S Eqd, k Peg, and V819 Her 



Star 



HJD- 
2400000.5 



<5RA 
(mas) 



<5Dec 
(mas) 



(/^as) 



^ maj 



<Pc 0"RA TDgc 

(deg) (fias) (^as) 



LDC Align 



Rate 
(Hz) 



Eclipse 



Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 

Equ 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
K Peg 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 
V819 Her 



53508.50939 
53550.41402 
53552.38996 
53571.33634 
53584.32135 
53586.30477 
53607.23078 
53613.22582 
53614.22329 
53637.20185 
53656.13045 
53874.48665 
53909.42854 
53957.31149 
53970.25612 
53971.25621 
53977.26493 
54003.21361 
54005.19180 
54028.10922 
54037.11296 
54230.50355 
54266.47236 
53494.50786 
53586.36471 
53637.29586 
53921.45172 
53963.35858 
53978.32708 
53995.24384 
54003.32618 
54008.31093 
54075.12667 
54265.39966 
53109.47951 
53110.48183 
53123.45772 
53130.44208 
53137.43044 
53144.42426 
53145.39541 
53168.33949 
53172.35221 
53173.33202 
53181.33391 
53182.33164 
53186.30448 
53187.30462 
53197.26851 
53198.24258 
53199.29186 
53208.25236 
53214.24077 
53215.23094 
53221.22209 
53228.20946 
53229.22073 
53233.18295 
53234.20151 
53235.21764 
53236.16733 
53249.16006 
53466.52265 
53481.50628 
53494.45305 
53585.24930 
53874.42827 
53956.21398 



-86.2658 
-94.5630 
-95.1309 
-99.1010 
-101.4805 
-102.0695 
-105.7105 
-107.1684 
-106.8636 
-111.2007 
-114.1207 
-133.8986 
-134.0082 
-133.6026 
-133.8368 
-133.8811 
-133.7524 
-133.0033 
-133.0818 
-132.3883 
-132.9388 
-116.0819 
-111.8997 
104.6067 
83.9658 
71.2323 
25.0960 
-9.1846 
-14.0804 
-17.0621 
-22.4829 
-21.6029 
-38.5846 
-84.8163 
49.6406 
48.0946 
49.1860 
48.4778 
48.3928 
47.7017 
48.3082 
47.0275 
47.3441 
47.1604 
46.4333 
46.5646 
45.6584 
46.1427 
46.1852 
46.2937 
44.0252 
46.4303 
45.6337 
45.6364 
46.2559 
45.0884 
45.2196 
45.0993 
44.8269 
45.2153 
44.4598 
44.2467 
31.4132 
30.4432 
29.5257 
23.9818 
1.2427 
-5.0810 



-117.6468 
-143.7433 
-144.7397 
-155.5409 
-162.5823 
-163.5727 
-174.9135 
-177.7125 
-178.5240 
-189.7814 
-198.7179 
-275.5257 
-283.8970 
-293.0248 
-294.8599 
-294.9923 
-295.9196 
-299.3443 
-299.4815 
-301.8627 
-302.3034 
-301.5397 
-297.2144 
43.3687 
55.9066 
62.9289 
85.9815 
98.3661 
99.4673 
101.1081 
101.1592 
102.0991 
107.2133 
119.9685 
-84.4966 
-84.1334 
-85.9318 
-86.4135 
-87.1396 
-87.6612 
-87.8013 
-89.7513 
-90.1337 
-90.3599 
-90.7857 
-91.0136 
-91.1213 
-91.2225 
-92.1522 
-92.2555 
-91.9446 
-92.4901 
-93.3429 
-93.5172 
-92.9726 
-94.3310 
-94.5160 
-94.8458 
-94.7637 
-94.9183 
-94.8862 
-95.8032 
-102.8881 
-103.3348 
-103.1671 
-102.6116 
-88.7776 
-81.9129 



25.7 

33.6 

10.0 

9.8 

10.4 

13.2 

3.1 

6.7 

5.9 

12.3 

16.3 

12.3 

6.1 

9.0 

5.7 

7.3 

5.3 

16.7 

11.6 

13.3 
11.0 
13.1 
13.2 
13.6 
2.2 
6.3 
43.1 
6.2 
2.7 
11.5 
10.9 
4.6 
7.0 
8.2 
7.3 
11.9 
18.1 
6.6 
14.0 
25.1 
13.7 
15.0 
6.3 
8.0 
7.5 
13.9 
18.2 
13.0 
4.8 
5.7 
24.7 
6.6 
5.5 
4.8 
8.8 
7.2 
6.4 
6.0 
7.6 
8.5 
4.7 
4.3 
9.9 
11.1 
18.6 
10.1 
6.8 
9.2 



1872.8 

248.4 

85.5 

68.1 

104.4 

831.5 

80.6 

243.8 

174.7 

357.4 

215.1 

424.2 

58.2 

405.5 

50.6 

62.3 

39.3 

206.4 

115.3 

406.9 
319.6 
508.2 
112.0 
645.9 
11.5 
53.5 
10007.6 
79.3 
17.5 
514.2 
849.5 
114.1 
183.8 
389.0 
282.5 
600.4 
507.8 
205.8 
280.2 
1039.6 
316.5 
339.9 
170.0 
77.4 
174.6 
333.2 
426.3 
441.5 
117.1 
54.6 
1645.4 
181.6 
125.9 
110.6 
342.2 
100.6 
80.0 
64.7 
37.8 
107.1 
78.1 
87.9 
204.3 
311.1 
251.9 
114.7 
99.5 
341.7 



150.49 
151.98 
150.13 
149.71 
151.49 
150.94 
148.50 
150.28 
150.30 
157.38 
154.15 
147.53 
152.62 
154.61 
153.56 
152.11 
157.55 
161.12 
157.21 
153.38 
159.14 
146.70 
158.12 
143.23 
166.51 
173.82 
160.71 
165.54 
162.83 
159.38 
3.09 
2.48 
3.06 
142.94 
158.77 
159.53 
162.47 
162.94 
164.34 
167.13 
161.59 
162.93 
168.29 
33.97 
169.71 
169.62 
166.80 
166.94 
164.87 
160.37 
171.42 
37.67 
169.45 
167.53 
38.91 
169.45 
172.84 
167.67 
172.74 
176.57 
166.59 
172.71 
163.01 
38.18 
163.98 
174.02 
168.07 
170.23 



1630.0 

219.8 

74.3 

59.0 

91.9 

726.8 

68.8 

211.8 

151.8 

330.0 

193.8 

357.9 

51.7 

366.4 

45.3 

55.1 

36.4 

195.4 

106.4 

363.9 
298.7 
424.8 
104.1 
517.5 
11.2 
53.2 
9445.7 
76.8 
16.8 
481.3 
848.2 
113.9 
183.5 
310.4 
263.4 
562.5 
484.3 
196.7 
269.9 
1013.5 
300.3 
325.0 
166.5 
64.4 
171.8 
327.8 
415.1 
430.1 
113.0 
51.5 
1627.0 
143.8 
123.8 
108.0 
266.3 
98.9 
79.4 
63.2 
37.5 
106.9 
76.0 
87.2 
195.4 
244.6 
242.2 
114.0 
97.4 
336.8 



922.7 

120.4 

43.5 

35.4 

50.7 

404.0 

42.2 

121.0 

86.7 

137.9 

94.9 

228.0 

27.3 

174.1 

23.1 

29.8 

15.8 

68.6 

45.9 

182.7 
114.3 
279.2 
43.5 
386.8 
3.4 
8.5 

3306.4 

20.7 

5.8 

181.5 

47.1 

6.7 

12.0 

234.5 

102.5 

210.3 

153.9 

60.7 

76.8 

232.9 

100.8 

100.8 

35.1 

43.8 

32.1 

61.6 

99.0 

100.6 

30.9 

19.1 

246.7 

111.1 

23.7 

24.3 

215.0 

19.7 

11.8 

15.0 

8.9 

10.6 

18.7 

12.0 

60.4 

192.5 

71.8 

15.6 

21.6 

58.7 



-0.99949 
-0.94899 
-0.96437 
-0.94762 
-0.97257 
-0.99930 
-0.99617 
-0.99794 
-0.99697 
-0.99531 
-0.98163 
-0.99794 
-0.96785 
-0.99836 
-0.96133 
-0.96095 
-0.93151 
-0.96621 
-0.96173 
-0.99668 
-0.99468 
-0.99843 
-0.94550 
-0.99904 
-0.75093 
-0.66476 
-0.99990 
-0.95159 
-0.87056 
-0.99770 
0.97271 
0.73308 
0.81528 
-0.99904 
-0.99707 
-0.99819 
-0.99240 
-0.99360 
-0.98202 
-0.99386 
-0.98964 
-0.98778 
-0.98309 
0.97548 
-0.97114 
-0.97328 
-0.98197 
-0.99110 
-0.98715 
-0.94836 
-0.99488 
0.99718 
-0.97194 
-0.97962 
0.99860 
-0.92813 
-0.83914 
-0.91188 
-0.51352 
-0.60015 
-0.96552 
-0.93121 
-0.98524 
0.99731 
-0.96316 
-0.76046 
-0.94665 
-0.98729 



263 
370 
1352 



1089 



774 
639 



4856 



1601 

1597 

570 

526 

264 



2646 



1127 
2387 
1682 
3320 
835 
1614 
632 
556 
933 
1315 
702 
5532 



1 
1 
1 
1 
1 
1 
1 

1 
1 
1 





1639 



615 
2244 
5863 
389 
1590 
2850 
1756 
3661 
2011 
1334 
1378 
2537 
1226 
897 
1673 
1409 
2560 
2904 
2795 
2014 
706 
1578 
5218 
5404 
946 
6558 
5251 
5723 
3998 
3180 
3905 
3303 
3701 
2094 
6684 
5428 
3031 
3301 
1355 
1479 
2358 
4028 1 















100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

100 

50 

50 

50 

50 

50 

50 

50 

50 
100 
50 
50 
100 
100 
100 
50 
100 
50 
100 
50 
50 
50 
50 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
50 



Note. — All quantities arc in the ICRS 2000.0 reference frame. The uncertainty values presented in these data have not been resealed. Column 1 is the star name. Columns 2-14 
are as columns 1-13 in Table^ Column 15 is if the measurement was taken during a subsystem eclipse, 1 otherwise {V819 Her only). 



TABLE 5 
Orbital parameters for, 5 Equ 





Value 


Uncertainty 


P (days) 


2084.03 


±0.10 


T (MHJD) 


53112.071 


±0.052 


e 


0.436851 


±0.000025 


a (mas) 


231.9650 


±0.0080 


V(),Licfc (kms-1) 


-15.40 


±0.11 


V(),DAO (kms-i) 


-15.875 


±0.080 


Vo,c (kms-i) 


-15.73 


±0.10 


Ml (Mq) 


1.192 


±0.012 


M2 (Mo) 


1.187 


±0.012 


Ml + M2 (M0) 


2.380 


±0.019 


M1/M2 


1.004 


±0.012 


i (deg) 


99.4083 


±0.0098 


u) (deg) 


7.735 


±0.013 


n (deg) 


23.362 


±0.012 


d (pc) 


18.379 


±0.048 


TT (mas) 


54.41 


±0.14 


Note. — All parameter uncertaint 




crcajscd by a factor of 


■Jx?^ = 1.09 (though the of 


the combined fit is artificial due to 


rescaling the un- 


certainties of the indiv 


dual data sets 


this reflects the 


degree to which the data sets agree 


with each other). 


The fit was repeated s 


everal times varying the set of 


non-dcgcncratc paramcl 


crs used in ordc 




tainty estimates for a n 




d quantities. The 



parameters {a, R = M1/M2} were replaced with the sets 
{M = Ml +M2,R} and {Mi, M2}. The parallax is a 
derived quantity. 



TABLE 6 

Orbital parameters for k Peg and V819 Her 







Peg 


V819 Her 


Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


Pab (days) 


4224.76 


±0.74 


2019.66 


±0.35 


Tab (MHJD) 


52401.52 


±0.96 


52627.5 


±1.3 




0.3140 


±0.0011 


0.67974 


±0.00066 


lAB (degrees) 


107.911 


±0.029 


56.40 


±0.13 


OJAB (degrees) 


124.666 


±0.064 


222.50 


±0.22 


Qas (degrees) 


289.037 


±0.021 


141.96 


±0.12 


PBaBb (days) 


5.9714971 


±1.3 X 10-6 


2.2296330 


±1.9 X 10-6 


TBaBb (MHJD) 


52402.22 


±0.10 


52627.17 


±0.29 


esoSi) 


0.0073 


±0.0013 


0.0041 


±0.0033 


iBaBb (degrees) 


125.7 


±5.1 


80.70 


±0.38 


i^BaBb (degrees) 


179.0 


±6.0 


227 


±47 


^BaBb (degrees) 


244.1 


±2.3 


131.1 


±4.1 


Vo,Keck (kms-1) 


-9.46 


±0.22 






Vi^Keck (kms-^day-i) 


-2.2 xlO" 


4 ±3.4 X 10-4 






V2,Keck (kms-1 day-2) 


6.8 X 10-8 


±2.4 X 10-6 






Vo,c (km s-i) 


-9.40 


±0.26 






Vo,Lick (km s-l) 


-8.37 


±0.26 






Vo,M/K (km s-l) 






-3.375 


±0.059 


Vo,DAO (km s-i) 






-3.385 


±0.065 


Vo,DDO (km s-l) 






-3.35 


±0.12 


Ma (Mq) 


1.533 


±0.050 


1.799 


±0.098 


MBa+Bb (Mq) 


2.472 


±0.078 


2.560 


±0.067 


MBb/MBa 


0.501 


±0.049 


0.742 


±0.012 


LiBb/LBa 


0.015 


±0.021 


0.280 


±0.037 


d (parsecs) 


34.57 


±0.21 


68.65 


±0.87 


*AB-BaBb (degrees) 


43.4 


±3.9 


26.3 


±1.5 


MBa (Mq) 


1.646 


±0.074 


1.469 


±0.040 


Msb (Mq) 


0.825 


±0.059 


1.090 


±0.030 


aAB (AU) 


8.122 


±0.063 


5.108 


±0.046 


O-BaBb (AU) 


0.08710 


±0.00091 


0.04569 


±0.00040 


7r (mas) 


28.93 


±0.18 


14.57 


±0.19 



Note. — Uncertainties for k Pegasi are the maximum of three uncertainties; 
the uncertainty from the combined fit that included PHASES-reweighted data, that 
including PHASES data with a 161/ias noise floor, and the difference in the fit values 
for the two models. The parameters are the average values between a fit including 
the reweighted uncertainties and one with the noise floor. 
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